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Abstract 

The photon energy spectrum in inclusive weak radiative B —* X s ~f decay is 
computed to order a^Po. This result is used to extract a value for the HQET 
parameter A from the average (1 — 2£^ 7 /ms), and a value of the parameter 
Ai from ((1 — 2E 1 /tub) 2 ) ■ An accurate measurement of (1 — 2E 1 /m,B) can 
determine the size of the nonperturbative contributions to the T(IS') mass 
which cannot be absorbed into the b quark pole mass. 
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Comparison of the measured weak radiative B — > X s ~/ decay rate with theory is an 
important test of the standard model. In contrast to the decay rate itself, the shape of 
the photon spectrum is not expected to be sensitive to new physics, but it can nevertheless 
provide important information. First of all, studying the photon spectrum is important for 
understanding how precisely the total rate can be predicted in the presence of an experimen- 
tal cut on the photon energy jj]], which is important for a model independent interpretation 
of the resulting decay rate. Secondly, moments of the photon spectrum may be used to 
measure the heavy quark effective theory (HQET) parameters which determine the quark 
pole mass and kinetic energy 0|§, much like the shape of the lepton energy Q or hadronic 
invariant mass spectrum in semileptonic B — > X c I v decay. The main purpose of this pa- 
per is to present the order a^j3 piece of the two-loop correction to the photon spectrum, and 
to study its implications. A calculation to this order is required for a meaningful comparison 
of the HQET parameters extracted from B — > X s 7 with those from other processes. 

To leading order in small weak mixing angles the effective Hamiltonian is 
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where Gf is the Fermi constant, Vij are elements of the Cabibbo-Kobayashi-Maskawa ma- 
trix, Ci(fi) are Wilson coefficients evaluated at a subtraction point fi, and 0j are the dimen- 
sion six operators 
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In Eq. (0), e is the electromagnetic coupling, g is the strong coupling, m;, is the b quark mass, 
is the electromagnetic field strength tensor, G* is the strong interaction field strength 
tensor, and T a is a color SU(3) generator. The sums over q include q — u,d, s, c, b and 
the subscripts L, R denote left and right handed fields. The Wilson coefficients have been 



calculated to next-to-leading order (NLO) Using a s (rriz) = 0.12, and the convention 

that the covariant derivative is = d^ + igA^T 0, + ieQA^ (where Q is the fermion's electric 
charge), the values we need are C^m^) = 1.13, C-j{m b ) = —0.306, Cs(m b ) = —0.168 ||. 

For the photon energy, E 1 , not too close to its maximal value, the photon spectrum 
dT/dE^ for weak radiative B decay has a perturbative expansion in the strong interaction 
fine structure constant a s . It is known at order a s and the main purpose of this letter is 
to present the order a 2 fto (so-called BLM H) contribution. It is well known that the part 
of the order a 2 piece proportional to the one-loop beta function, j3 = 11 — 2n//3 usually 
provides a reliable estimate of the full order a 2 piece. This part of the order a 2 s contribution 



is straightforward to compute using the method of Smith and Voloshin [|T(| . 

Using the dimensionless variable[|, x b = 2E^/m b , the photon energy spectrum in B 
X s ~f takes the form 

a s (m b ) . , . , / a s (m b y 
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where 
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is the contribution of the tree level matrix element of O-j to the B — + X s ^f decay rate, and 

Ci(m b )Cj(m b ) 



A p (x b ) = J2a^(x b ) 



(5) 
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The sums over i,j in Eq. (|5|) give the contributions of the various operators in Eq. (0) to 
the photon energy spectrum. 

It is important to note that since the coefficients in H e g are known only to NLO accuracy, 
the BLM calculation of the 0\ — 0$ contribution to the photon spectrum is only meaningful 
away from the endpoint. At the endpoint, order a 2 contributions to the matrix elements are 



1 Later we will introduce a dimensionless photon energy variable normalized by the B meson mass, 
xb = 2E 1 /niB- 
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the same order as the unknown NNLO running [where a s hi(mw/mb) is counted as 0(1)]. 
Neglecting the small contribution to A from 0\ — 0§ discussed in the next paragraph, at 
least one gluon must be in the final state to populate the spectrum for xj, < 1, so it is 
consistent to combine the a 2 matrix elements with the NLO Wilson coefficients. (Strictly 
speaking, we should for consistency only use the (3q part of the NLO running of the operators 
with the BLM calculation, but for simplicity we will use the full NLO result. The difference 
between these two approaches is small.) Thus powers of a s in Eq. @ and elsewhere reflect 
the perturbation expansion of the matrix elements only, and not of the Wilson coefficients. 

At zeroth order in the strong coupling, the spectrum for xj, < 1 arises from matrix 
elements of the four-quark operators 0\ — O e in Eq. (||). Of these 0\ and 2 include two 
charm quarks in the final state, and therefore they contribute to the photon spectrum only 
for lower values of x& than what we consider in this paper. These contributions are divergent 
in perturbation theory, and the divergence can be absorbed into the definition of the quark to 
photon fragmentation function, D 9 ~* 7 (x), which depends on an infrared scale A. D q ^(x) is 



calculable in the leading logarithmic approximation JTT| , |r2| . There is some data on D 



X 



however, the experimental errors are still quite large |13|| . This fragmentation contribution to 
the coefficients clq(x) vanishes as xt, — > 1, and it is small in the region of large Xb, 0.65 < x&, 
which we consider in this paper. 

A very important B decay background to the B — ► photon spectrum is from non- 
leptonic h — > cud and b — > uud decays, where a massless quark in the final state radiates 
a photon. Such backgrounds due to the operators (ci , y tJl bi J ){di J n ffj,ui) and ( y ULl^bi J )(di J 'y^UL) 
are shown in Fig. 1 (using IV^/V^I = 0.1). We used the Duke-Owens parameterization of 
the fragmentation function ||14|| , setting A = 1.3 GeV and Q 2 = mf. (This value of A is 



motivated by a fit to the ALEPH data |I3"|.) The uncertainty of this result is sizable, since 
the A-dependence is large and may not be large enough to justify keeping only the lead- 
ing logarithms. Close to maximal Xb the resummed fragmentation function may predict too 
large a suppression of the photon spectrum, since the lightest exclusive final states dominate 
there. The background from b cud (b — > uud) is more than 50% of the 77 contribution 
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FIG. 1. B decay background to the photon spectrum due to the operators {cLl^bL){dLl '^ul) 
(solid curve) and {uL")^bL){dL^ '^ul) (dashed curve). 

to (1 /I 'o)dT I 'dxb below Xb ~ 0.75 (xb ~ O.65).0 Therefore, we will concentrate on the region 
Xb > 0.65; to measure the B — > X s ~f photon spectrum at lower values of Xb would not only 
require excluding final states with charm with very good efficiency, but also demanding a 
strange quark in the final state. Note that for B — > Xd'j, the fragmentation contribution 
from b — > uud is larger than the short distance piece unless Xb is very close to 1. 

Neglecting the strange quark mass, a 88 is also divergent in perturbation theory. This 
divergence can also be absorbed into the definition of fragmentation functions. In the leading 



logarithmic approximation |T5 
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where D s ^ J (x) and D g ^ 1 {x) are the strange quark to photon and gluon to photon fragmen- 
tation functions, which have large uncertainties. In the region Xb > 0.65, the a 88 contribution 
to the photon spectrum (l/T )dT/dxb is less than 0.01. Given the uncertainty in a 88 , and 
its small magnitude, it does not appear useful to calculate a 88 . 



2 Note that these backgrounds are steeply falling functions of Xb, and are indeed negligible in the 
present CLEO region of E~ > 2.1 GeV. The tree level contribution of the operators O3 — Oq in 



Eq. (||) to the photon spectrum is about a fifth of the b — > uud background. 



Experimentally, because of backgrounds, only B — > X s 7 photons with large energies can 
be detected. The present experimental cut is E 1 > 2.1 GeV at CLEO 0, which corresponds 
to xi, > 0.875 with m& = 4.8 GeV. In the large Xj, region the most important contribution to 
the sum in Eq. (|5[) come from the 77 term, with moderate corrections from the 22, 78, and 
27 terms. The other contributions (88, 28, and the ones involving 0\ and O3 — Oq) are very 
small, and will be neglected in this paper. 

Simple analytic expressions for a[ 7 and a] 8 are available, 

(2x 2 - 3s - 6)2 + 2(x 2 - 3) ln(l - x) 



a\\x) 
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Neglecting the small Aq term in Eq. (|3[), we can calculate the shape of the photon spectrum 
away from x — 1 to order a 2 /3o accuracy knowing the effective Hamiltonian to order a s 
(NLO) only. At order a 2 /3o, we find that a 1 ^ and a 18 are given by 



a 7 2 7 (x) 
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where L 2 {z) = — Jq dtln(l — t)/t is the dilogarithm. The strange quark mass is neglected 
throughout this paper; it only enters the final results quadratically, as m 2 /[m 2 (l — Xf,)]. 

The functions of a 22 and a 27 are known in the literature P^JT/| , and we agree analytically 
with those results. The order a 2 j3o contributions, a 2 ? and computed numerically. 

We find it most useful to present simple approximations to these functions 
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-0.0842 + 0.3333a; - 0.2005a; 2 + 0.0227a; 3 
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and 




-0.1064 + 0.4950a; - 0.4361a; 2 + 0.0373a; 3 




-0.0156 + 0.0463a; + 0.3467a; 2 - 0.3045a; 3 + 0.0027 ln(l 

I 771 1 4\ 
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These approximations are accurate to within 1% in the region x b > 0.6 for m c /m b = 1.4/4.8. 
The 27 contribution is very sensitive to m c /m b . Changing m c /m b from 1.4/4.8 to 1.2/4.6 or 
1.6/5.0 modifies af and a 27 dramatically. The 22 contribution only changes in the previously 
mentioned range of m c /m b by ±(20 — 25)%. The 22 contribution is also accurate to within 
1% when m c /m b changes by ±0.03. However, the 27 contribution is only accurate at the 
20% level when m c /m b changes in this range. Note that the perturbation series in a s is 
particularly badly behaved for the 27 contribution. Roughly 2/3 of the 22 contribution is 
from absorptive parts corresponding to real intermediate states. 

The coefficients are determined by calculating the order a 2 n/ piece and making the 
identification, — 2n//3 — > (3q. There is a subtlety in applying this method to weak radiative 
B decay. There is a contribution of order a° s rif from the tree level b — > s^qq matrix elements 
of O3 — Oq, coming from Feynman diagrams where the photon couples to the bottom or 
strange quarks. It is not associated with a term of order a°/9 . To avoid adding an analogous 
spurious order a 2 /9 contribution to af and a 22 , only diagrams where the photon couples to 
the charm quark were included in the calculation of the matrix element of O2. 

Part of the B — > X s 7 matrix element of O2 is not adequately calculated in perturbation 
theory. It corresponds to the process B — > J/ipX s followed by the decay J ftp — > 7 + 
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FIG. 2. The sum of the 77, 22, 78, and 27 contributions to (l/To)dT/dxb at order a s (thick 
dashed curve) and ai?/3o (thick solid curve). The thin curves show the 77 contribution only. The 
scale is the same as in Fig. 1. 

(light hadrons). There will be large corrections to the part of the charm quark loop where the 
cc are almost on-shell and have the same velocity. In this region there are large "Coulombic 
QCD corrections" that produce the J/ip state. However, cutting this small part of the cc 
phase space out of our calculation of the matrix element of O2 has a negligible effect. Hence, 
at the order of perturbation theory to which we are working, calculating the cc loop while 
removing J/^'s from the data would be a consistent approximation. 

The sum of the 77, 22, 78, and 27 contributions is plotted in Fig. 2 in the region 0.65 < 
Xb < 0.9 (using a s {m b ) = 0.22 and /3 = 25/3). For very large x, other effects that we have 
not calculated become important. There are both nonperturbative and perturbative terms 
that are singular as x — > 1. They sum into a shape function that modifies the spectrum 



in this region [|ll|. Unfortunately, at the present time, it is not possible to make a model 
independent estimate of these effects. Therefore, we do not plot the perturbation theory 
predictions for x b > 0.9. In the plotted region, the 22, 78, and 27 terms make a moderate 
correction to the dominant 77 contribution to (l/T )dT/dx, which is shown in Fig. 2 with 
the thin curves. 

The b quark mass can be eliminated in favor of the B meson mass by a change of variables 

to 
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xb = 2E 1 /m B ■ 



(13) 



Using m b = m B — A + (Ai + 3X 2 )/(2m b ) + . . ., the photon spectrum becomes 

dT 



dx 
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1 + + 

m B 



dT 
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x b =x B (l+A/m B +...) 



(14) 



For xb within a region of order Aqcd/^b of unity (its maximal value) nonperturbative 
effects are very important. However, for integrals of xb over a large enough range these 
nonperturbative effects are small. 

An important integral of this type is 

dT 

(15) 



[ dx B (l-x B ) 
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The parameter 5 = 1 — 2E™ in /tub has to satisfy 6 > Aqcd/^b; otherwise nonperturbative 
effects are not under control. It is straightforward to show that 
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where 



f 1 1 

l-x b ) = / dar 6 (1 - ar 6 ) — 



i dr 

dxh 



(17) 



Note that all terms but the first one in Eq. (J16|) have perturbative expansions which begin 
at order a s . The ellipses denote contributions of order (Aqcd/^b) 3 , as(AQCD/^s) 2 , and 
a 2 s terms not enhanced by @o, but it does not contain contributions of order (Aqcd/^s) 2 
or additional termsQ of order a s (Aq CD /m B ). Terms in the operator product expansion 



3 There are actually additional contributions formally of order o s (Aqcd/"ib) coming from the 
expansion of m c /m b in the 22 and 27 terms. Although the 27 term is very sensitive to the value of 
mc/mh, this A-dependence is negligible for (1 — xb)\ Xb >i-8- 
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FIG. 3. The sum of the 77, 22, 78, and 27 contributions to (1 — Xb)\ Xb >i-S at order a s (thick 
dashed curve) and a 2 /3o (thick solid curve). The thin curves show the 77 contribution only. 

proportional to Xi^/ml en t er precisely in the form so that they are absorbed in mg in 
Eq. ([16]) ||. There are also nonperturbative corrections suppressed by (Aq CD /m c ) 2 instead 
of (AQ CD /m b ) 2 |19| . These do not contribute to Eq. ([161) . 

Using our results, (1 — %b) \ Xb >i~s in Eq. (|T7| ) is known to order a 2 /3 . Writing 



;i - xb) 



x b >l—& 



IT \ 7T / 



(18) 



B p have decompositions analogous to Eq. 



C^m^C^rrib) 



C 7 2 (m 6 ) 



(19) 



Neglecting B (5), Eqs. (0) and @ yield for the dominant 77 contribution 

6^(5) = — [-95 3 + 145 2 + 725 - 54 + 12(5 2 - 35 - 6) In 6 
54 L 



(20) 



&-(5) 



- 3695 4 + 1165 3 + 18005 2 - 38525 



2592 

+ 408tt 2 + 125(95 3 + 345 2 - 1025 + 66) In 5 

- 2165(5 2 - 35 - 6) In 2 5 - 144(5 3 - 35 2 - 65 + 17)L 2 (1 - 5) 



(21) 



Our prediction for (1 — Xb)\ Xb >is is shown in Fig. 3 as a function of 5, both at order a s 
and a 2 /?o- The bad behavior of the perturbation expansion would improve somewhat by 
evaluating the strong coupling at a smaller scale than mj, such as rribVS, the maximal 
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available invariant mass of the hadronic final state. This bad behavior may also be related 



to the renormalon ambiguity in A. 

A determination of A is straightforward using Eq. ([16]). The left hand side is directly 
measurable, while (1 — Xb)\x b >i-s and (l/T Q )dT / dxb\ Xb >is in the second and third terms 
on the right hand side can be read off from Figs. 3 and 2, respectively. Using the CLEO 
data in the region E 1 > 2.1 GeV [l|], we obtain the central values A a 2^ ~ 270 MeV and 
A as ~ 390 MeV. We have indicated the order kept in the perturbation expansion to determine 
A, since a value of A extracted from data can only be used consistently in predictions valid 
to the same order in a s . These values are consistent with the ones obtained from a fit to 



the B — > X c lv lepton spectrum ||], and from the CLEO fit pTJ to the B — > X c £u hadron 
mass distribution 0. 

At the present time this extraction of A has large uncertainties. The potentially most 
serious one is from both nonperturbative and perturbative terms that are singular as x — > 1 
and sum into a shape function that modifies the spectrum near the endpoint. A model 
independent determination of these effects is not available at the present time, however, 
it may be possible to address this issue using lattice QCD |22|. For sufficiently large S 
these effects are not important. It has been estimated that they may be significant even if 



the cut on the photon energy is lowered to around E 1 = 2 GeV p3| , |24| , but this is based 
on phenomenological models. We have implicitly neglected these effects throughout our 
analysis. The validity of this can be tested experimentally by checking whether the value of 
A extracted from Eq. ([16]) is independent of 5 in some range. This would also improve our 
confidence that the total decay rate in the region i B > 1 — 5 can be predicted in perturbative 
QCD without model dependence. 

The value of A at order a s has a sizable scale dependence: lowering the scale such that 
a s changes from 0.22 to 0.3 reduces the value of A as by about 40 MeV. At order ai?A) this 
scale dependence is much smaller. Uncertainties due to the unknown order (Aqcd/^b) 3 
terms in the OPE [TM] are largely uncorrelated to those in the analyses of the lepton energy 



or hadron mass spectra in B — > X c £ v |25| . The effect of the boost from the B rest frame 
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0.35 



FIG. 4. Prediction for (1 — xb)\ Xb >i-8 i n the upsilon expansion at order e (thick dashed curve) 
and (€ 2 )blm (thick solid curve). The thin curves show the 77 contribution only. 



into the T(AS) is small for (1 — xb)\ Xb >i- 



The upsilon expansion [T2^] yields parameter free predictions for (1 — xb)\ Xb >i~s in terms 
of the Y(IS') meson mass. The analog of Eq. fllED is 



x b >1 — <5 



fflT 

2m 



B 



1 + O.Olle + 0.019(e 2 ) BL M - (1 - x b ) 



Xb>(2m B /m?)(l-6) 



(22) 



where e = 1 denotes the order in the upsilon expansion. For E 1 > 2.1 GeV this relation 
gives 0.111, whereas the central value from the CLEO data is around 0.093.0 In Fig. 4 we 



plot the prediction for (1 — xb)\ Xb >i~s as a function of S, both at order e and (6 2 )blm- The 
perturbation expansion is much better behaved than the one shown in Fig. 3. The most 
important uncertainty in this approach is the size of nonperturbative contributions to the 
T(1S) mass other than those which can be absorbed into the b quark mass. These have 
been neglected in Eq. (p2|). If the nonperturbative contribution to the T(IS') mass, A T , were 
known, it could be included by replacing mx by mx — At- For example, Ax = +300 MeV 



4 It is interesting to note that including the CLEO data point in the 1.9 GeV < E~> < 2.1 GeV 
bin, the experimental central value of (1 — xb) over the region E y > 1.9 GeV is 0.117, whereas the 
upsilon expansion predicts 0.120. 
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FIG. 5. The sum of the 77, 22, 78, and 27 contributions to ((1 — Xb) )\ Xb >\-5 at order a s (thick 
dashed curve) and a^Po (thick solid curve). The thin curves show the 77 contribution only. 



increases (1 — xb) by 21%, so measuring (1 — xb) with such accuracy will have important 
implications for the physics of quarkonia as well as for B physics. 

The variance of the photon energy distribution can be used to determine Ai The 
analog of Eq. (|i~6"D in this case is 



x b >1-8 



1 -XB, 



1 2 



x B >l— 6 



3m B ~^ 3 ^ ' 



(l-— )((1-^) 2 ) 
V m B J 



tub 1 o dxb 



Xh=l—S 



x b >l—8 

+ ..., (23) 



where (3 ~ 0.064 is the magnitude of the velocity of the B meson in the T(4S) rest 
frame, and only the leading /3-dependence has been kept. The ellipses denote terms of 
order (Aqcd/^b) 3 , «s(Aqcd/^_b) 2 , and a 2 s terms not enhanced by (3q. Our prediction for 



((1 — Xb) 2 )\ Xb >i-5 is shown in Fig. 5. Note that unlike the case of (1 — xb)\x b >is, the effect 
of the boost is very important in Eq. fl2"3"D- Using the CLEO data in the region E 1 > 2.1 GeV, 
we obtain the central value Ai ~ —0.1 GeV 2 , with large experimental errors. The uncertainty 



in this value of Ai due to A is small. Nonperturbative effects from the cut on E 1 p4| , and 
the unknown higher order contributions to Eq. (|23|) are expected to have a larger impact on 
the determination of Ai than the corresponding effects have on the determination of A from 
Eq. (H). 
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In summary, we calculated order a 2 s j3o corrections to the shape of the photon energy 
spectrum in weak radiative B — > Xgj decay. The dominant 77 contribution is given by 
simple analytic formulae in Eqs. (|7|) and (|]). The other terms relevant in the region Xb > 0.65 
are the 22 and 27 contributions given in Eqs. ( [TT]) and (|T2"D, and the 78 term given in Eqs. (||) 
and (]TD|) . The HQET parameter A can be extracted from the average (1 — 2E 1 /mg) using 
Eq. and it can also be used to test whether the nonperturbative contribution to the 

Upsilon mass is small. The CLEO data in the region E 1 > 2.1 GeV implies the central 
values A as ~ 390 MeV and A Q 2^, ~ 270 MeV at order ot s and a^Po, respectively. Possible 
contributions to the total decay rate from physics beyond the standard model are unlikely to 
affect this determination of A. In the future, checking the 5-independence of the extracted 
value of A, and comparing the experimental and theoretical shapes of the photon spectrum 
for Xb < 0.9 can provide a check that nonperturbative effects and backgrounds are under 
control. This would also improve our confidence that the total decay rate in the region 
xb > 1 — 5 can be predicted model independently, and used to search for signatures of new 
physics with better sensitivity. 
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